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Abstract 

The systems engineering description of a wideband communications channel is provided 
which is based upon the fundamental propagation aspects of the problem. In particular, the well 
known time variant description of a channel is formulated from the basic multiple scattering 
processes that occur in a random propagation medium. Such a connection is required if optimal 
processing methods are to be applied to mitigate the deleterious random fading and 
‘multipathing’ of the channel. An example is given which demonstrates how the effective 
bandwidth of the channel is diminished due to atmospheric propagation impairments. 


Introduction 

From the communications engineering standpoint, the intervening propagation channel 
between a transmitter and receiver is treated as a separate system-element whereby only the input 
and output characteristics of the channel are defined; the intricate structure and mechanism 
within the system-element communications channel is not identified nor considered (see, for 
example, [ Rummler , 1979]). The associated atmospheric transfer function is traditionally 
empirically detennined. Such an approach has found applicability in the design and specification 
of a communications system that must operate in the deleterious channel mechanism of the 
earth’s atmosphere (see, for example, [ Proakis , 1989]). However, in the case of wideband 
transmission through planetary atmospheres which can display a host of different propagation 
mechanisms, one needs to add a descriptive level to the channel which reflects these 
mechanisms. Thus, one must go back to fundamentals and ‘re-derive’ the system engineering 
description of a communications channel from basic propagation principles. From such a 
detailed description can one effectively specify and design optimal processing algorithms to 
mitigate the prevailing multipath and fading of the propagation channel. This goal forms the 
subject matter of this work. 


Fundamental Propagation Aspects of the Problem 

The relevant starting point is, of, course, the Maxwell equations; with a spatially and 
temporally dependent random permittivity e = s(r,t), one has for the electric field E = E(r,t ) 
and magnetic field B = B(r,t ) of the modulated wave field 
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which is a stochastic vector wave equation due to the random source term on the right side. 

Equation (4) is a general relation describing the evolution of the vector electric field 
propagating through the random permittivity field £ . The term on the right side of this equation 
describes the depolarization of the electric field as it propagates thus making eq. (4) a vector 
equation for the three components of the electric field [Manning, 1993, p. 8]. Although a general 
solution can be obtained for this equation, this circumstance complicates the subsequent analysis 
using this approach and a simplification is possible whereby the depolarization term is neglected. 
The necessary conditions under which this is valid have been established and can be found 
elsewhere [ Manning , 1993, appendix A]. So as not to deter form the present analysis, this 
depolarization term will be neglected in what is to follow. 

Thus, with these considerations prevailing, the neglect to the depolarization term allows 
eq. (4) to decouple into three scalar equations, one for each component to the E field. Therefore, 
one is now concerned with the stochastic scalar wave equation 

vi£ -^F (f£)=0 <5) 
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for one of the components E of the electric field. At this point, it is important to remember that 
the permittivity field e(r,t) is a random function of position r and time t so this equation does 
not have any straightforward solution and must be treated only with statistical methods; with the 
E field being a random field, due to the presence of the e function, the E field can only be 
described in terms of its statistics, i.e., its mean, variance, etc., in space and time. Departing from 
the usual “narrow-band” treatment (see appendix A) of eq. (5) by which it is reduced to the 
Helmholtz equation, the field E will be taken to contain a range of spectral frequency 
components due to modulation of the propagating wave field; thus, eq. (5) must be dealt with in 
its entirety. 

In what is to follow, the function e will be written as a sum of a deterministic part 
£ 0 = e 0 (r,t) = 1 and a random part e = e(r,t) with the associated conditions |e| « 1 and (e) = 0. 
(Here, the notation (y) indicates the ensemble average of a random quantity y.) Making this 
assignment in eq. (5) and simplifying gives 



Using the temporal Fourier Transform 


E(r,co) = J E(r,t)exp(icot)dt 


( 7 ) 


and applying it to eq. (6) gives 
( 1 \ 2 1)2 


1 Y T d 2 ( 1 T d~ 

V 2 £’(r,m)- — J — jE(r,t)exp(icot)dt = — J — j[e{r,t)E{r,t)]exp(i(Ot)dt (8) 

\c J _^dt yc J ^dt 

The second term on the left side of this equation becomes, after two integrations by parts with 
applicable boundary conditions at t — » , 


f ] V 


1 f T 

— J — jE{r,t)exp(i(Ot)dt = —k 2 E{r,(6) 
dt 


Similarly, the member on the right side becomes 
( 1 V 2)2 


— I J ^-j[e(r,t)E(r,t)\exp{icot)dt = -k 2 J e{r,t)E{r,t)exp(i(Ot)dt 
V c / -oodt 


Equation (8) can then be written in the temporal frequency domain as 


( 9 ) 


( 10 ) 


V 2 E(r,co)+ k 2 E(r,co) = -k 2 J e(r,t)E(r,t)exp(icot)dt 


( 11 ) 
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The exact relation given by eq. (11) can be simplified by using the quasi-stationary 
approximation (see appendix A) whereby, relative to the inherent temporal variations of the 
functions E(f,t ) and exp (icot) that contribute to the integrand of eq. (11), the temporal variations 
of e(r,t) are ‘slow’ (see eq. (A. 3) of appendix A). Thus, one can approximately factor £(?,/) out 
of the integrand on the right side of eq. (11) and obtain the simpler relation 

V 2 E(r,co) + k 2 E(r,co ) = -k 2 e{r,t)E{r,(6) (12) 

which is isomorphic to that used in investigations of the propagation of monochromatic 
radiation. 

Any of the well-known methods of analysis in stochastic wave propagation theory can now 
be used to study the ramifications of this relation [Manning, 1993], i.e., the Rytov method, the 
parabolic equation method, or many of the methods developed for quantum field theory. In order 
to maintain generality in attempting to relate the rigorous methods of wave propagation theory 
with the ad hoc methods employed in communications theory, the widely applicable quantum 
theoretic methods will be used. In particular, use will be made of a factorization technique which 
will concatenate within a single parameter all orders of multiple scattering of the wave field 
within the random medium, resulting in what is known as the Lippmann-Schwinger equation; the 
single parameter is known in quantum field theory as the T -matrix but, for purposes here will be 
called the random scattering operator T . With this in hand, it will be shown that the resulting 
relations can indeed be put into the form of time-variant systems theory which is so widely used 
in multipath propagation but, more importantly, giving explicit meaning to the quantities of the 
time-variant systems model. 

To this end, one writes the solution of eq. (12) as (see appendix B) 

E(r,eo) = E 0 (r-R,ct)j+ J G 0 (r,r')q(r',co)d 3 r' (13) 

where E 0 (r - R, ft>) is the initial field distribution due to a source at R of frequency ft), 

q{r\(o) = -k 2 e{r\t)E{r',co) (14) 

is the scattered field source term and the Green function G 0 (r,r') is the solution of the related 
differential equation 


VG 0 (r,r') + k 2 G 0 (r,r') = 8{r - r') 


(15) 


which is 


G 0 (r,r') = G 0 (\r-r' |) 


1 cxp(/f 

r - 

' F 1) 

An 

r - 

- r' 



(16) 


Thus, remembering that k = co/c which allows eq. (16) to be stated in terms of ft), the solution 
for eq. (12) becomes 
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e(r',t)E (r',0))d 3 r 


(17) 


r 


k 2 T CXP| 

E(7,co) = E 0 (r-R,co) + — J 


co ^ 
i — \r — r\ 


\r — r 


Hence, even though the propagation of modulated wave fields is considered in this 
development, due the fact that the spectral content of the wave field is significantly smaller than 
the frequency of the carrier wave, and due to the quasi-stationary approximation, the form of the 
solution of the randomly scattered field is identical to that of the monochromatic case. Although 
eq. (17) will be developed further in what is to follow, a brief digression will be made 
concerning the meaning of the content of eq. (17) over that of the monochromatic case which 
will form the foundation of the analysis of digital transmission. 

Returning to the time-domain using the Fourier transform relationship inverse to that of 
eq. (7), viz., 


E(r,t) 



(18) 


and applying this to eq. (17) yields 


E(r,,)=E 0 {?-Rj) + E± J J 


= E 0 (r-R,t} + 


4k 2k 

kr_ 
4k 


exp| 

(.CO 
l — 

\r — r‘ 

] - icot\ 
£ 


— — / 

r — r 



— J exp| -ia\ t- 
2k 


\r — r\ 


€(?, t)E{r',(6)cEr'd(0 

i!£7) dV 


E{r\co)clco 


V 


\r — r 


= E ^-Rj) + E] E 


r,t- 


\r~r'\)£(r\t) d3/ 


) 


r — r 


(19) 


In the Fraunhofer diffraction region, in which the receiving point is sufficiently removed form 
the scattering position, one can expand the quantity | r - r'|, viz., 

\r — r \ ~ r — n ■ r -\ — 

where n = rjr is the unit vector directed from the center of the scattering volume to the receiver 
point thus allowing eq. (19) to be written as 


E(r,t) 





( 20 ) 


where V is the volume comprising the scattering region. Finally, iterating once the solution of 
eq. (20) yields the Born approximation for the single scattering events, i.e., 
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( 21 ) 


E(r,t) 


E «(r-R.t)+j^:fe(7',t) E 



Multiple Scattering Within the Channel 

The form of the solution in eq. (21) shows how the resulting total field E(r,t ) is the sum of a 

direct, unscattered field, i.e., E 0 (r-R,tj, and a scattered component which is a function of the 

time delay t d = r/c - h ■ r'/c. Such multipath phenomena results in a temporal spread induced 
into the communications channel; see, e.g., [ Proakis , 1989]. In addition, the perturbation 
introduced by the random permittivity e(r,t) will give rise to a fading channel. 

In fact, if one were to continually iterate the solution given by eq. (20), a multiple scattering 
series is obtained whereby the n th term in the series describes n multiple scattering events. Such 
multiple scattering gives rise to a total received signal which is a sum of signals that each 
traverse separate paths through the propagation medium resulting in multipath interference. This 
phenomena serves to degrade the available bandwidth of the communications channel. This will 
be shown in a much more explicit form in what is to follow. 

Returning now to eq. (17) with the goal of obtaining a relationship that includes all orders of 
multiple scattering, one commences with the straightforward iterative solution of eq. (17). At the 
outset, however, it will be beneficial to simplify the notation somewhat. Defining the ’propagator’ 
function 


P(r,r',CQ ) 


. CO r 

2 exp i — \r — r\ 
k L c 

An \r — r'\ 


( 22 ) 


allows eq. (17) to be rewritten as 


E(r,a > ) 


E 0 (r-R,coj+ J P(r,r , ,^o)e(r , ,^)£ , (r , )<r/ 3 r , 


(23) 


As is shown in appendix C, one can obtain an iterative solution to eq. (23) which, upon 
factorization, can be written as 


E(r,a > ) 




J P(rJ-\a)T(r\r",t,a)E 0 (r"-R,(d 



/ 

r 


(24) 


where the factor 


T(f',r",t,(o ) = e(r',n 8(7'- 7") + P{f',r'',(6)e(7",t)+ J P(r',7 m ,co)£(r m ,t)P(7' ir J-'',co)£(r'',t)di" + 


(25) 
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is the scattering operator which convolves all possible multiple scattering events that occur 
between the scattering inhomogeneities e . 

Equation (24) is of the form of the well-known Lippmann-Schwinger equation of quantum 
electrodynamic scattering theory [Lippman and Schwinger, 1950], [ Goldberger and Watson, 
1964], As mentioned above, the utility of this relation rests in the explicit connection of the 
initial, unperturbed field E 0 with that of the total scattered random field E through the use of the 
random scattering operator T convolving all the scattering mechanisms within the random media 
through which the wave propagates. 

Equations (24) and (25) can also be written in a convenient but less explicit operator 
notation, i.e., 


E = E 0 + PTE 0 , T = e(l + Pe + PePe + • • •) (26) 

where the T in the expression for the scattering operator is the integral identity S(r — r'). The 
operator representation for T suggests a further factorization allowing it to be written as the 
recursive operator relationship 


T = e( 1 + PT) 


(27) 


This is indeed the case as can be established by using the same factorization method as employed 
to derive eq. (24) which gives in coordinate representation 


T(r,r',t,0)) = e(r,t ) S(r — r') + J P(r,r",a)T(r\r'',t,6))d 3 r" 


(28) 


The Time-Variant Systems Description of the Channel 

Equation (24) is of a form that easily lends itself to the establishment of a time-variant 
systems description of propagation through a random medium. To this end, one lets the initial 
field distribution be given as a temporally modulated wave 

E 0 (r ) = A(r)m(t)exp(-ico 0 t) (29) 

where A[r ) is the initial spatial distribution of the wave field, m(t) is the modulating signal, and 
( o 0 is the carrier frequency. Using the Fourier transform relationship of eq. (7) gives 

E 0 (r,co) = A(r)M(co-CO 0 ) (30) 


where 


M(co-co 0 )= J m{t)exp[i(co-co 0 )t]dt 


(31) 
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is the frequency spectrum of the modulation. Substituting eq. (30) into eq. (24) yields 


E(r,O)) = M((O-(O 0 )v{r-R,CO,tj 


( 32 ) 


where 


V(r -R,C0,t^j = A(r - R^+ J J P[r,r\(o)T(r',r /, ,Qj,t)A(r /, —Rjd 3 r , d 3 r (33) 


Inverse transforming eq. (32) back to the temporal coordinate gives 


jj M(v)v{r - R.v + m 0 , / j c x p (-/ w ) r/ v c x p (-/ m 0 / ) 


Defining the envelope e(r,t) of the field by the relation E(r,t) = e(r,t)exp(-ico 0 t) gives 

e {r,t)=( 


M(v)v(r-R,v + ® 0 ,t)exp(-zw)r/v 


Using the relation 


(34) 


(35) 


M(v ) = J m(t)exp(z w)<r/t 


(36) 


in eq. (35) and rearranging terms allows one to finally write for the envelope of the modulated 
propagating wave field 


e(r,t) = J m{tyi(r — R,t — (37) 

where the impulse response of the atmospheric communications channel is given by 

h(r - R,t - t',t) = — J vlr - R,v + CQ 0 ,t)exp[-iv(t - t^dv (38) 

From this impulse response function, one can define correlation as well as power spectral 
density functions from which one can define several characteristics of the fading multipath 
communications channel [ Proakis , 1989]. For example, one can form expressions for the delay 
power spectrum, the spaced-frequency, spaced-time correlation function, as well as the Doppler 
power spectrum of the channel. However, the most direct demonstration of the effects of the 
multipath and fading on the modulation attendant on the channel is through the response function 
for the average power. 


NAS A/TM— 2004-212975 



One can obtain the response function for the average power as follows. It is important to note 
that the quantity v(r — R,v + co 0 ,tj is a random function since the scattering operator which 
enters into its definition (eq. (33)) is a random function of the permittivity e(r,t) given by eq. 
(25). Thus, any variable or function related to it will also be a random quantity. Thus, one can 
only refer to the measurable aspects of these quantities in terms of their statistical descriptions 
such as the average and higher order statistical moments. For this reason, one can only calculate 
the average power from the use of eq. (35), viz, 

I(r,t) = (e(r,t)e*(r,tfj 

= J J M(v)M\v^(v(7-lv + o) 0 j)v\7-lv'+co 0 ,t))- 

•exp [~i(v - v^tjdvdv' (39) 


Defining the two-frequency mutual coherence function (MCF) by 

r(r - R, v + c o 0 , V + (O 0 , ij = (v{r-R,v + co 0 , ijV* [r - R,v'+ co 0 , tj 


(40) 


as well as the center and difference frequencies, respectively, 

v + V 


C = 


v d = v-v 


(41) 


allows eq. (39) to be written as 


7 M= J J M V c +U?^y* V c ~~ v d jr(r -^,v c + (O 0 ,v d ,t^jexp[-iv d t]dv c dv d 


(42) 


where the MCF is taken also to be a function of the center and difference frequencies of v + c o 0 
and v' + co 0 . 

This general relation can be simplified if the spectral densities of the modulation are such 
that they fall-off as a function of v c for values of v c « (O 0 thus essentially yielding a negligible 
contribution to the integral at values of v c ~ co 0 . In this instance, one can write 

r(r-/?,v c + O) 0 ,V d ,tj ~T (r-R,co 0 ,v d ,ij, i.e., the MCF becomes approximately independent of 
v c and, thus, eq. (42) takes the form 






r-R,a> 0 ,v d ,t) exp [-i v d t ] cl v d 


(43) 


where 
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(44) 


1 M 




Ml 



is the input modulation power spectrum and the MCF r( r - R,co 0 ,v d ,tj acts as a time-variant 
transfer function. 

Proceeding further in the simplification of this result, one can apply a quasi-steady state 
assumption based on the fact that the channel characteristics, ultimately connected to the time 
varying permittivity s(r,t), slowly evolve over the time scale of digital transmission along the 
link. In this event, the MCF can be made independent of time allowing one to write 
T (r -R,co 0 ,v d ,tj ~r(T-R,O) 0 ,v d y Equation (43) then becomes 


7 M=k- 



L )(' 


-R,(D 0 ,v d 


)exp [~iv d t]dv d 


(45) 


Fourier transforming this relation with respect to co yields 


I(r,co)= r M 



(46) 


This shows how, in the case for a time-independent channel in which the MCF is a function of 
only the difference frequency, the modulation power spectrum (m/2) is modified by the MCF 


A Simple Time-Invariant Example 

An illustrative example will now be given which shows how, due to the deleterious multipath 
effects concatenated in the time-independent MCF r(f,m 0 ,m/2), the effective bandwidth of the 
time-invariant communications channel is degraded. Consider the propagation of a Gaussian 
pulse given by 


3 


m(t) = expl 


o J 


From eq. (36), the associated spectrum of the pulse is 


M(v ) = VflT 0 ex p 


( ^ \ 


hL v 2 


(47) 


(48) 


The bandwidth Av associated with this pulse spectrum can be operationally defined by 
exp(-l) = exp^-r o 2 (Av)~/4) ; thus 
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Av = 


(49) 


_2 

T 

J o 


The input signal power spectrum corresponding to eq. (48) is, from eq. (44), 



(50) 


Now suppose, for ease of calculation, that the atmospheric MCF of the channel is that of a plane 
wave (i.e., independent of spatial position) also of Gaussian form, i.e., 



(51) 


where the coefficient f3((Q 0 ) is a function of the propagation parameters that characterize the 
propagation path, viz., number and type of aerosols, amount of atmospheric turbulence, etc. 
Using eqs. (50) and (51) in eq. (46) and performing the inverse Fourier Transform gives 



(52) 


The situation with f5((O 0 ) = 0 corresponds to the ‘vacuum’ case in which the effects of the 
intervening atmosphere are absent. In this instance, eq. (52) gives for the vacuum case 




-exp 


2r 

T 2 

V 7 o ) 


(53) 


Comparison of these two expressions yields the fact that the effective bandwidth of the 
communications channel is decreased from that given in eq. (49) to that given by 


Av = 



(54) 


Summary 

The usual time-variant systems theory approach to the design and analysis of 
communications systems that must operate through a random channel, such as that of the 
atmosphere, has been derived from first principles using multiple scattering theory. In particular, 
the atmospheric transfer function of the communications channel is directly related to the 
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multiple scattering equations of stochastic propagation theory. From this, the two-frequency 
mutual coherence function can be defined by considering the average power transmitted within a 
modulated wave field. It was then shown how, in the simple time-invariant case, this two- 
frequency coherence function modifies the modulation power spectrum within the random 
channel and results in a degraded channel bandwidth. 
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Appendix A 


At the outset of all ‘narrow-band’ analyses of the stochastic Maxwell equations used to 
describe statistical propagation situations [Manning, 1993, pp. 2-7, and references therein], one 
assumes a wave field E = E(r,t ) of single temporal frequency ft) and decomposes it using the 
relation 


E(r,t ) = E(rj)exp(-ift)/) 


allowing eq. (5) to be written as 


V 2 E + k 2 eE = 


2ik <9(eE) 1 <9 2 (eE) 

c dt c 2 dt 2 


(A.l) 


(A. 2) 


where k = co/c is the free space wave number. In this expression, the wave field E = E(r,f) is still 
a random function of time but is devoid of the temporal oscillatory component. The only 
temporal variation enters in through the random permittivity function e(r,t). One writes this 
function as a sum e(r,t) = £ 0 (r,t) + e(r,t) of a deterministic part £ 0 (r,t) and a random part 
e(r,t), and for simplicity, lets £ 0 (r,t) = 1 and assumes |e(r,t)| « 1. Using this subsequent 

relation in eq. (A. 2) and analyzing the partial temporal derivatives that appear shows that they 
can be deemed negligible so long as the condition 


d£ 1 T <9E \ 

dt Ef dt j 


(A. 3) 


prevails. In this instance, eq. (A. 2) simplifies to 

V 2 E + k 2 £ E = —k£ E 


(A. 4) 


i.e., a stochastic Helmholtz wave equation which is only parametric in time.. The condition of eq. 
(A. 3), which states that the level of the temporal variation of the permittivity field is much 
smaller than that which is characteristic of the electric field of the propagation wave, is thus 
known as the quasi-static approximation. The structure of eq. (A. 4) is such that the random 
function e, which can be statistically specified depending upon the propagation scenario, 
multiplies the function being sought, i.e., the stochastic field E. This circumstance is the source 
of difficulty in obtaining solutions to eq. (A. 4) and fonns the subject of the many investigations 
on the subject in the past 40 years. The approach taken in the present work assumes an arbitrary 
range for the frequency content of the wave field necessitating the retention of the temporal 
derivatives in eq. (5). It is important to remember that in eq. (5), both the functions £ and E are 
functions of time and thus the indicated operations contain several derivative terms which, in the 
present case, must each be examined for their level of contribution. 
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Appendix B 

To easily show this development, one writes, for notational clarity, eq. (12) in operator 
notation, viz., 

LE(r,co) = q(r,co), L = 'V 2 +k 2 (B.l) 

and assumes a solution of this differential equation to be of the form 

E(r,<x>) = E 0 ^-R,co}+ \ G 0 (r,r']q(r',(o)d 2 / (B.l) 

where E 0 (r - R,a)j is the initial (unscattered) field distribution and the kernel G 0 (r,r') of the 
integrand must be determined from the structure of the problem given in eq. (B.l). Thus, 
applying the operator L to eq. (B.2) and using eq. (B.l) gives 

LE(r,co ) = LE 0 [r - R,0)j + L J G 0 (r,r^(r',co)d 2 r' 

= J LG 0 (r ,ryj(r',cd)d 3 / 

= q(r,(o) (B.3) 

where, by definition, one has LE 0 (r - R. to] = 0 , and for the last two lines of this development to 
hold, one must necessarily have 

LG 0 (r,r ') = S(r - r') (B.4) 

which is eq. (15). The solution of eq. (B.4) can be obtained using Fourier Transform methods. 
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Appendix C 

Iterating successive expressions of eq. (23) gives (suppressing the t and (0 dependence for 
simplicity) 

E(?) = E 0 (7-R ] )+ } P(r,ry(^E 0 (r'-R)d 3 r + 

+ } } P(Ery{ry{r,ry(r^E 0 (r-R)d 3 r'd 3 r''+ 

+ } } ]p(rjy{ry(r,ry£{r^P{r\r^£(r)E 0 ^ m -R)d 3 r'd 3 r , 'd 3 r m +--- (C.l) 

The basic idea here is to factor out appropriate quantities within each of the integrals within each 
term on the right side so as to isolate those factors that exclusively describe the multiple 
scattering process. To this end, one can immediately factor an integral and a propagator function 
allowing eq. (C. 1) to be written 

E(r) = E 0 (r-H) + J P (?,?'){§(?') E 0 (?' - R) + 

+ } P(r,r , )£(?'')E 0 (?''-R)d 3 r''+ 

+ J } P(f(f^£(f^P(f',f^£(f^£ 0 (f w -^)dVVr w +-- ]</V (C.2) 

One now needs to factor out the initial field distribution E 0 [r - ft) ; it is thus required to isolate 

each of the factors E 0 {r - R^j in each member of eq. (C.2) so they can be factored from the 

remaining portions of the integrands; the latter will then form a series in multiple scattering of 
the wave field within the random medium. Making use of the fact that the dummy variables in 
each of the integrands can be interchanged with one another as well as, from eq. (22), 
P(r,r ') = P(r',r), the fourth member of the right side of eq. (C.2) can be rewritten as 

} } P(r,ry(r^P(r\r m )£(r^E Q (r-RYr r d 3 r m = 

= } } P(f(f^£(f*)P(f w ,f^£(f')£ 0 (f'-^)^V^ 3 r w 

where the transcription r" <-> r'" has been applied. Thus, the former factor E 0 (r"' - has been 
changed to E 0 \r" — Rj without changing the result of the double integral; this process can, in 
principle, be carried out for all remaining terms that appear in the series. It is very important to 
note that this process can only be employed to convert all C w (f (,i) - /?) factors to E 0 {r" - R^. 

This is so because if each such tenn were converted to E 0 (r' - R^j, i.e., with a single prime on the 
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r coordinate, the third member on the right of eq. (C.2) would not yield to the process due to the 
lack of enough dummy variables. Thus, attention must now be turned to the second member on 
the right side to convert the E 0 (r' - R^j factor to one involving the coordinate r" , i.e., 

E 0 \r‘ " - Rj. One can simply use the prescription E 0 [r' - R^j - S(r' - r")E 0 [r" - Rj noting that 
the integral over the r" coordinate will return the proper function. With these considerations in 
mind, E 0 [r" - R^j factors can be taken out of each term along with the integral over r " , and eq. 
(C.2) finally becomes 

E(?) = E 0 {r-R ] j+ J J P(r,ry , (r\r^E 0 (r'-R)d 3 r'd 3 r (C.3) 

where 

T(r',r') = £{?') 5(7-7") + P(7,7")e(7")+ } P(7, 7 , y(7')P(7' , J-y(?y 3 7'+--- (C.4) 
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